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We study the thermal conductivity within the E\ g and E2 U models for superconduc- 
tivity in UPts and compare the theoretical results for electronic heat transport with 
recently measured results reported by Pussier, Ellman and Taillefer. The existing 
data down to T/T c ps 0.1 provides convincing evidence for the presence of both line 
and point nodes in the gap, but the data can be accounted for either by an E\ g or E' iu 
order parameter. We discuss the features of the pairing symmetry, Fermi surface, 
and excitation spectrum that are reflected in the thermal conductivity at very low 
I/") temperatures. Significant differences between the E\ g and E' iu models are predicted 

to develop at excitation energies below the bandwidth of the impurity -induced An- 
dreev bound states. The zero -temperature limit of the c axis thermal conductivity , 
limy^o k c /T, is universal for the E' iu model, but non-universal for the E\ g model. 
C Thus, impurity concentration studies at very low temperatures should differentiate 

between the nodal structures of the E' iu and E\ g models. 
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Introduction 

Recent experiments of the thermal conductivity at low temperatures in UPt3 
have narrowed the range of possible symmetry states for the order parameter. 
Among several proposed order parameter models the two-dimensional orbital rep- 
resentations coupled to a symmetry breaking field appear to be very promising 
candidates. 1-4 For UPt3, which has a hexagonal crystal structure (Dgh), phase dia- 
gram studies 5,6 and transport measurements 7-9 lead to either an even-parity, spin 
singlet Ei 3 or an odd-parity, spin triplet F<2 U pairing state. In this paper we com- 
pare theoretical results for the electronic thermal conductivity of these two models 
with recent experimental data on single crystals of UPt3. 10, 11 We use the quasi- 
classical theory of superconductivity to calculate the thermal conductivity with the 



formulism derived in Ref. [12]. Our study focuses on the low temperature limit 
where the heat transport is dominated by electronic quasiparticles scattering off 
impurities. 

In order to develop a theoretical model that accurately predicts the low temper- 
ature behavior of clean anisotropic superconductors we must determine the basis 
functions for the order parameter. These functions are the eigenfunctions of the 
linearized gap equation and are therefore determined by the microscopic pairing 
interaction. Since the pairing mechanism is uncertain for the heavy fermion super- 
conductors we are forced to adopt a phenomenological approach in which calcula- 
tions of the transport coefficients, based on models for the pairing basis functions, 
are compared with experiment. Fortunately, crystal symmetry provides significant 
constraints on the basis functions. The eigenfunctions of the linearized gap equation 
form irreducible representations of the point group. Thus, the ground state order 
parameter and excitation gap often exhibit nodal lines or points on the Fermi sur- 
face, which reflect the symmetry class of the order parameter. 13 This is the case for 
the 2D representations that have been proposed as models for the ground state of 
UPt3. The basis functions for the pairing states that are often used in calculations 
are spherical harmonic functions or tight-binding functions with the transformation 
properties of the appropriate representation. However, it is worth emphasizing that 
the basis functions are not uniquely defined by symmetry. Thus, quantitative anal- 
ysis is expected to require more detailed knowledge of the pairing functions than is 
contained in the simplest models for the basis functions. Our analysis focuses on the 
low temperature limit for two reasons: (i) the transport coefficients exhibit univer- 
sal features in the limit T —> 0, which reflect the symmetry of the pairing state, and 
(ii) the low temperature transport coefficients are determined by lower dimensional 
regions of the Fermi surface near the nodes of the gap function. Thus, the complete 
Fermi surface and the detailed shape of the excitation gap are not required in order 
to accurately calculate the transport properties at low temperatures, T<T C . 

In Ref. [12] it was shown that certain eigenvalues of the thermal conductivity 
tensor are universal (i.e., independent of the scattering rate) below a characteris- 
tic temperature set by the bandwidth of impurity-induced Andreev bound states. 
Whether or not a universal limit develops at low temperatures depends sensitively 
on the nodal structure of the excitation gap. Both the Ei 3 and F<2 U models for the 
order parameter lead to an excitation gap with a nodal line in the basal plane and 
point nodes along the ±c directions. 14 The difference in the excitation spectrum for 
these two states is that the gap opens linearly for small angles away from the c di- 
rection for the Ei 3 order parameter, but quadratically for the F<2 U order parameter. 
This difference is reflected in limy^o k c /T, which is universal for the F<2 U model, 
but non-universal for the Ei 3 model. A detailed description of these asymptotic 
limits is given in Ref. [12]. 

Recently Fledderjohan and Hirschfeld 15 examined the temperature dependence 
of the anisotropy ratio, for the Ei 3 and F<2 U models and concluded that the 

data of Lussier et al. 10 favored an Ei 3 state with point nodes vanishing linearly 
near the poles. The authors qualified this conclusion because their calculations 
assumed a spherical Fermi surface and lowest order spherical harmonics for the 
pairing basis functions. Very recently, Norman and Hirschfeld 16 carried out numer- 
ical calculations of the electronic thermal conductivity using a Fermi surface for 
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UPt3 constructed from local density approximation (LDA) calculations, and sev- 
eral models for the pairing basis functions (e.g., tight-binding functions and Fermi- 
surface harmonics). These authors found considerable improvement between theory 
and experiment for the absolute magnitudes of Kj_ y/T compared to the results in 
Ref. [15] for both Ei 3 and F<2 U models. They reported a good fit between theory 
and experiment for the Ei 3 state, but they did not find an adequate fit in the F<2 U 
case. However, we note that the calculations reported in Ref. [16] do not describe 
the data of Lussier et al. 10 accurately in the low-temperature region, T < 0.3T C , for 
either pairing symmetry; the data of Lussier et al. 10 for k/T does not exhibit the 
saturation shown in the theoretical curves of Ref. [16], which assume a scattering 
rate of order To = 0.1 T c . This same observation was made by Lussier et al. 10 and 
led them to conclude that the low temperature data for kj /T is more compatible 
with To ~ 0.01 T c and a pairing symmetry of E2 U . However, all of these conclu- 
sions are based on calculations with model basis functions that do not accurately 
reproduce the low energy excitation spectrum. 

The E-Rep Models 

In our comparison of theory with experiment we pursue a different approach 
than Norman and Hirschfeld. 16 Instead of examining the effects of the multi-sheet 
(LDA) Fermi surface on the heat current, we model the excitation spectrum by 
an excitation gap that opens at line (ab plane) and point (c axis) nodes on the 
Fermi surface, and by the Fermi surface properties in the vicinity of the nodes 
(i.e., the Fermi velocities near the nodes). The low-temperature behavior of the 
transport coefficients probes 'lower-dimensional' regions of the Fermi surface, where 
the excitation gap vanishes, and is less sensitive to the overall geometry of the Fermi 
surface. It is therefore admissible to parametrize the nodal regions of the gap, and 
attempt to fit the nodal parameters in order to achieve accurate low temperature 
limits for the thermal conductivity along the principal axes of the crystal. The 
advantage of this approach is that we can quantitatively determine the phase space 
contributing to the low temperature transport coefficients and then examine in more 
detail the effects of impurity scattering and order parameter symmetry on the heat 
current, without having to know the overall shape of the Fermi surface or basis 
functions. 

We consider two models for the order parameter in the low-temperature phase 
of UPt3 which belong to different representations of the Deh symmetry group: 1 (1) 
the even-parity Ei 3 (hybrid-I) pairing state, and (2) the odd-parity F<2 U (hybrid-II) 

state, 

(Eig) A(pj) = A (T)y Elg ($), (1) 
(E 2u ) A(pf) = A (T)y E2 Jp f )c, (2) 

where c is the quantization axis along which the pairs have zero spin projection. 2 

1 Since we are interested in the low-temperature phase of UPt3 we neglect the small effect of the 
symmetry breaking field. 

2 In this analysis spin has no effect on the transport coefficients beyond the connection between 
the direction of the spin quantization axis and the nodal structure of the Fj2 U basis functions. 
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The corresponding basis functions are 

yE lg (ff ) = Vo Pfz (Pfx + ipfy ) (1 + a 2 pj z + a 4 pf z ) , (3) 
yE 2u (Pf ) = 2/o Pfz {Pfx + ipfyf (1 + a 2 pf z + a 4 pf z ) , (4) 

where 2/0 is chosen to normalize the basis functions such that Ao is the maximal 
gap. The coefficients a 2 and a 4 determine the variation of the gap near the nodes. 
In the vicinity of the equatorial line (pf z = 0), |A(#)| ~ /fii ne Ao|0|, with //i; ne = 
?/o, while near the poles (p fx = p fy = 0), |A(6»)| ~ Anoint A |6|" with fi pomt = 
(1 + a2 + 04) 2/0, and n = 1 for Ei 3 and n = 2 for E2 U . 3 This parametrization 
allows us to adjust independently the opening of the gap (slope or curvature) at 
the line and point nodes. For a 2 = a 4 = the standard spherical harmonic basis 
functions are recovered. The crucial difference between the Ei 3 and F, 2u states lies 
in the opening of the gap with angle at the polar point nodes, as shown in Fig. 1. 
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Fig. 1. The excitation gap (exaggerated) at the Fermi surface (sphere) for Ei 3 (a) 
and (b) order parameters showing the line node in the basal plane and the 
linear and quadratic point nodes. Panel (c) shows the normalized excitation gap at 
T = for the same two states as a function of the polar angle 0. The parameters 
correspond to those that fit the low temperature thermal conductivity data shown 
in Figs. 2 and 3. 

3 The parameters fi\' me and /ipoint d enne the slope or curvature of the gap near a nodal line or 
point in a spherical coordinate system that is obtained by mapping an ellipsoidal Fermi surface 
onto a sphere. 
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The transport properties of unconventional superconductors are strongly influ- 
enced by scattering from impurities. One of the more striking effects is the appear- 
ance of a band of low energy excitations deep in the superconducting state. 17-19 
This occurs for non-magnetic impurities in unconventional superconductors when 
the Fermi surface average of the order parameter vanishes, (A(pf)} = 0, i.e., the 
order parameter changes sign around the Fermi surface. Impurity bound states 
develop from the combined effects of impurity scattering and Andreev scattering. 
In an unconventional superconductor with a line of nodes a finite concentration of 
impurities leads to a band of Andreev bound states with a bandwidth, 7 < Ao, be- 
low which the density of states is almost constant and non-zero at zero energy. 20,21 
This novel metallic band, deep in the superconducting phase, exhibits universal 
transport properties at very low temperatures. 12, 22-25 

In Ref. [12] it was shown that the principal components of the thermal con- 
ductivity tensor can be expressed in terms of an effective transport scattering time 
that incorporates all of the coherence effects of superconductivity for T —> 0, 

v 2 - 

KmKi(T) = ~~ Js T Ti (i = a,b,c). (5) 

Here v 2 i is the Fermi-surface average of [vf,i(pf)] 2 , Js = f 71 " 2 ^^/ is the normal- 
state Sommerfeld coefficient, and the effective transport scattering time is defined 

by 



v. 



2 



>J,i J [A(p f ) 2 + J 2 



(6) 



where 7 is the bandwidth of Andreev bound states. For heat flow in the basal 
plane one has r a = r b = ta, where ta = 37i/(4 / uii ne Ao(0)), which is universal, 
i.e., independent of the scattering cross section or impurity concentration. The 
leading temperature corrections to K b /T are given by a Sommerfeld expansion for 

T<7<7 1 c , 

Kb (T) ~ a b T + p b T 3 . (7) 

In the strong scattering limit the parameters a b and f3 b are related to the microscopic 
model parameters by 

V fb n I 7tT 2 R 

a b ~^-y S T A , (i b a b ~ — — — . (8) 
3 60 7 2 

For heat flow along the c-axis the thermal conductivity again has a Sommerfeld 
expansion at low temperatures, k c = a c T + /3 C T 3 . However, the c-axis coefficients 
depend sensitively on the winding number of the order parameter around the point 
nodes, 12 i.e., n = 1 for Ei 3 and n = 2 for E'ju, 

n „ v h . n ., f 2 7/^ pomt Ao(0), (E lg ) 



3 7S 2// pomt Ao(0) 1 1, (E 2u ) W 

7* 2 k 2 B f 5 , (E lg ) 



M«< * W*t ^ (Su) • (10) 
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Note that a c = limj^o k c /T is universal for the F<2 U state, but depends on the 
impurity concentration for the Ei 3 state. In the unitarity scattering limit 

■y 2 oc Ti{, thus, the coefficient of the T 3 term scales with the impurity concentration 
as f3 c ~ l/rii for the F<2 U state and f3 c ~ l/ v /nj for the Ei 3 state. 



Numerical Results 

We have carried out numerical calculations of the thermal conductivity, based 
on the formulism developed in Ref. [12]. These calculations coincide with the ana- 
lytic results for the very low temperature limit described above. A principal input 
to the numerical calculations is the total quasiparticle scattering rate, T(T), which 
is the sum of an elastic (impurity) and an inelastic (electron-electron) contribution 
(Matthiessen's rule). Estimates of the impurity scattering rate can be obtained 
from the residual resistivity or normal state thermal conductivity. The normal state 
scattering rate obtained from the thermal conductivity, T(T) = ^v 2 b js T/ kj(T), 
is found to be nearly isotropic (s-wave scattering), and dominated by electron- 
electron scattering above the superconducting transition for clean single crystals of 
UPt3. 10 In the temperature range, 0.5 K < T < IK, the thermal conductivity is well 
described by the total scattering rate, T(T) = 2t (t) = T 2 /T 2 ) ; 10 for the 

particular crystal of Ref . [10] the elastic and inelastic contributions are nearly equal 
at T c . The inelastic rate decreases rapidly at low temperature, even faster than T 2 
because of the opening of a gap over most of the Fermi surface. However, below the 
superconducting transition we model the electron-electron scattering self-energy in 
the superconducting state by the same phenomenological, temperature dependent 
T(T) as in the normal state. Thus, we neglect the reduction in T(T) due to the onset 
of superconductivity. At T = 0.1T C the inelastic rate is certainly less than 1% of the 
total scattering rate. Since our primary interest is in the low temperature behavior 
of the heat transport the error in the inelastic rate is insignificant. However, we do 
not expect as accurate a description of the experimental data in the temperature 
regime, T c /2 < T < T c , where inelastic scattering is not negligible. 

In Figs. 2 and 3 we plot theoretical results for [tti(T) T c ]/[Ki(T c ) T] for the 
Ei 3 (Fig. 2) and F<2 U (Fig. 3) models. Impurity scattering is assumed to be in 
the unitarity limit and the slope parameter of the line node is //i; ne = 2. The 
experimental data for UPt3 shown in Figs. 2 and 3 was obtained from Ref. [11] and 
is normalized by the normal-state value kj(T c )/T c k, 20mW/(K 2 cm) at T c ~ 0.5 K. 
The data is in excellent agreement with the theoretical curves for either symmetry 
class in the low temperature range, 0.1T C < T < 0.4T C . The theoretical curves 
saturate for T < T* ~ j to the low temperature limits described by Eqs. (7)- 
(10). The experimental data shows no evidence of saturation down to the lowest 
temperature; thus, T* < 0.1T C . In addition to this constraint, the universal limit of 
the thermal conductivity in a superconductor with an equatorial line node implies 
that the ratio 

lim^^ ^) (11) 
T^o Kb (T c )T 2 W i„eA (0) v ' 

scales linearly with the scattering rate, T(T C ), and is independent of the scattering 
phase shift, So- The dependence on T(T C ) in Eq. (11) comes from the normal-state 



6 



value, Kb(T c ) ~ 1/T(T C ), and the parameter Hiine = (1/ Ao)\dA(6) / d6\ no d e describes 
the opening of the gap at the line node. 



— r(T c )=0.03jrT t 
- r(T c )=0.02jrT t 




T/T c 

Fig. 2. The normalized thermal conductivity of an even-parity Ei 3 -state in the 
unitarity limit for heat flow along the b and c axes, and for two different scattering 
rates of T(T C ) = 2To. The slope parameters at the linear point nodes and at the 
line node were fixed to \i\ = 2/3 and hi = 2.0. The experimental data (symbols) 
were taken from Lussier et al. 11 and normalized to the value at T c ~ 0.5 K. The 
inset shows the normalized anisotropy ratio (k c / kj)/ (k c / kj)t c - 



The value Hiine = 2 is consistent with the constraints imposed by the inter- 
cept of [ki,(T)T c ]/[ki,(T c )T] at T = and the crossover temperature, T*, i.e., 
r /(jUi m eAo) < | • 0.08 and T* ~ 0.2 ^Tofiy me A < 0.1 T c . The value of T is ad- 
justed to give the best fit to the low-T region of the basal plane thermal conductivity. 
For To = 0.0l7rT c o ~ 0.03 T c the intercept becomes, [ki,(T)/T]/[ki,(T c )/T c ]t=o ~ 
0.02. This impurity scattering rate, that fits the low-temperature part of the ther- 
mal conductivity in the superconducting state, is smaller than that estimated from 
the normal-state transport and dHvA data, To ~ 0.1 T c — 0.2T C , 26 by at least a 
factor of 3. This discrepancy is considerably smaller than the discrepancy reported 
by several other authors. 11, 10, 15 Finally, the parameters defining the excitation gap 
near the point nodes are adjusted to fit k c (T)/T at T <C T c . These values are 
Hi =2/3 for the Ei 3 model and ^2 = 4 for the F<2 U model. 
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Fig. 3. The normalized thermal conductivity of an odd-parity E2 U -state in 
the unitarity limit for heat flow along the b and c axes, with a scattering rate of 
To = 0.0l7rT c o, and several choices for the curvature parameter at the point nodes, 
H2 = 4.0, 4.5, 5.0. The slope of the line node was fixed to nn = 2.0. The inset shows 
the normalized anisotropy ratio (k c / ki,)/(k c / kj)t c - 



Although the existing low temperature data for the thermal conductivity on 
UPt3 is equally well described by either an Ei 3 or F<2 U model, the ultra low temper- 
ature region offers a possibility of distinguishing these two models. The anisotropy 
ratio of the thermal conductivities K c j K b as a function of temperature is shown in 
the insets of Figs. 2 and 3. The zero temperature limit of this ratio for the F<2 U 
(hybrid-II) state is a universal value, independent of impurity scattering, 

Kc (T)/ Kb (T) m (E } (12) 

T-0 K c (T c )/ Kb (T c ) ^2 ^ >' { ' 

while this ratio in the case of the Ei 3 (hybrid-I) state depends on the impurity 
concentration and scattering strength through the impurity bandwidth 7, 

Kc (T)/ Kb (T) _ W7 
t™o Kc (T c )/ Kb (T c ) ~ ^Ao(O) (tjls) - (U) 
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In the unitarity limit 7 ~ ^/ro/Uii ne Ao(0). The slope and curvature parameters 
describe the opening of the order parameter at the line node (jii, /^u) and at the 
point nodes (jii,^)- Thus, impurity concentration studies of the zero-temperature 
limit (i.e., T < T*) of the anisotropy ratio should distinguish between these two 
models for the ground state. 4 

The analysis has so far assumed unitarity scattering by impurities. It has been 
argued theoretically that impurity scattering in heavy fermion metals is generally 
expected to be in the unitarity limit. 27,28 Fig. 4 shows the influence of the scattering 
phase shift 60 on the thermal conductivity at low temperatures. It is clear that UPt3 
is in the strong scattering limit with a scattering phase shift 60 ^ 80°. For strong 
scattering the crossover from the power law region to the universal limit occurs just 
below T* ~ 0.1T C . However, the crossover to the universal limit occurs at much 
lower temperatures for weak scattering, and is barely visible for the scattering phase 
shifts less than 60 = 60°. 




Fig. 4. The normalized thermal conductivity of the E2 U -state for heat flow in 
the basal plane (along the 6-axis), a scattering rate T(T C ) = 0.027rT c o, and various 
scattering phase shifts So- The slope and curvature parameters are nn = 2.0 and 
H2 = 4.0. For weak scattering, 60 ^ 60°, the crossover region to the universal 
behavior becomes very small, and barely observable in the figure. The experimental 

4 The same suggestion is made by Norman and Hirschfeld. 16 
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data (circles) are from Lussier et al. 11 

The crossover to the universal limit is not seen in the currently available ex- 
perimental data of Lussier et al. 11 The observation of saturation of K{/T below T* 
would be an important test of theory. 11, 12 Recent measurements of heat transport 
in UPt3 by Suderow et al. 29 and Huxley et al. 30 extend down to T w 0.03 T c . These 
measurements are in agreement with a crossover temperature T* 0.07 T c . These 
authors also report an asymptotic value for the basal plane thermal conductivity of 
[n a (T)/T]T^o ~ 0.4 mW/(K 2 cm), which is in remarkably good agreement with our 
theoretical value, [kj/T]t^o ~ 0.02 • 20mW/(K 2 cm), extracted from the numerical 
fit to the data of Lussier et al. 10 shown in Figs. 2 and 3. A similar estimate for the 
universal value can be obtained with Eq. (8) from normal-state data. For UPt3 it 
is reported that j s ~ 10 kJ/(K 2 m 3 ) 30 and v s>b « v StC /y/2& « 3.3 km/s. 31 ' 30 Thus, 
we estimate for the heat current in the basal plane [ki(T)/T]t-*o ~ 1 mW/(K 2 cm), 
which is of the same order as our previous estimate 12 and is comparable to the value 
of 0.4 mW/(K 2 cm) reported by Suderow et al. 29 

Commentary 

At low temperatures heat transport is dominated by elastic scattering of quasi- 
particles in low dimensional regions of the Fermi surface near nodes of the order 
parameter. As a result, the low energy spectrum can be determined accurately, 
without complete knowledge of the Fermi surface. Excellent agreement with exper- 
iment at low temperatures is obtained for both pairing states, Ei 3 and E2 U . The 
differences in the excitation spectrum for these states are predicted to be observable 
at very low temperatures, T < T* ~ 0.1 T c , when the impurity-induced Andreev 
bound states determine the heat transport. Measurements of the thermal conduc- 
tivity ratio, k c /ki,, at very low temperatures for various impurity concentrations 
should differentiate between the two models; the F<2 U model is predicted to have a 
universal ratio, while this ratio for the Ei 3 state is predicted to depend strongly 
on the concentration of impurities. At high temperatures, T ~ T c , our calculations 
show deviations from the experimental data for both models, although the basal 
plane heat transport is in slightly better agreement for the Ei 3 model with a scat- 
tering rate of T(T C ) ps 0.06 T c . Accurate calculations of the heat transport in this 
temperature region require knowledge of the excitation spectrum over the entire 
Fermi surface, as well as a self-consistent treatment of inelastic scattering. 
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